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Abstract. We compute the one loop anomalous dimension of the gauge invariant dimension two 

operator "fuf J d^x (A^^)^, where U is an element of the gauge group, by exploiting Zwanziger's 
expansion of the operator in terms of gauge invariant non-local n leg operators. The computation 
is performed in an arbitrary linear covariant gauge and the cancellation of the gauge parameter 
in the final anomalous dimension is demonstrated explicitly. The result is equivalent to the one 
loop anomalous dimension of the local dimension two operator (^5^)^ in the Landau gauge. 
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The last decade has witnessed an intense interest in the area of dynamical gluon mass gen- 
eration related to the understanding of the low energy properties of Quantum Chromodynamics 
(QCD). See, for instance, [DElElilElElITlElElIiniin]. One motivation for this has been the 
numerical evidence for the apparent condensation of a dimension two operator to explain the 
deviation of an effective coupling constant from the expected perturbation theory prediction, 
[121 El El E] • To reconcile the difference one can fit the data more accurately with a l/Q^ 
correction requiring a dimension two operator to balance the dimensionality of the momentum 
Q. Ordinarily one would expect a dimension four operator correction as the leading power cor- 
rection due to the gauge invariant operator (GJ^^)^ based on the field strength GJJ^ where a is 
the adjoint colour index. That a dimension two operator appears to emerge in this analysis, 
[121 [T3l [T^ I15j . is not inconsistent with a variety of other observations made over a period of 
years. Indeed [161 1171 ^^8] noted that the perturbative vacuum is unstable and the condensation 
of a dimension two operator would be energetically favourable. Earlier studies of potential gluon 
mass operators included a Coulomb gauge analysis of a dimension two operator, [19], as well as 
Cornwall's construction of a massive gauge invariant QCD Lagrangian which supports vortex 
solutions, [201 121j . However, one main theoretical objection to such dimension two operators 
is that the obvious naive choice, ^(A^)^, where is the gluon field, is clearly gauge variant 
and therefore not suitable for condensing in quantities involving gauge invariant objects. As the 
effective coupling constant of [HI [131 El [15] is gauge dependent, there is therefore no immediate 
reason why such a gauge variant operator condensate cannot be the explanation of the 0(1/(5^) 
power correction. However, phenomenological analyses of gauge invariant quantities, [1], appear 
to require a gluon mass, albeit tachyonic, to fit experimental data. Motivated by Curci and 
Ferrari's work from the 1970's, [22], there has been a recent re-examination of ^(^^^ and its 
BRST invariant extension since their model possesses a classical gluon mass term. The main 
drawback of the Curci-Ferrari model for being a possible Lagrangian of a massive vector boson 
is that the BRST operator in not nilpotent and hence unitarity is absent, [231 124| . Instead the 
hope would be that the quantum condensation of the operator would circumvent the unitarity 
objection. Though this is still an open question. 

In this respect one calculation of note was the Landau gauge analysis of [1[ for the condensa- 
tion of \{A'^)'^ in Yang- Mills theory using the local composite operator method. This was later 
extended to QCD for massless quarks in [25j. The premise of investigation of [4] rests in the ob- 
servation that one can in fact have a gauge invariant dimension two operator in QCD which can 
condense. Indeed the resultant phenomenology of the particular operator in question has been 

elaborated on in [26]. This operator is given by m J d^x where U is an element of the 

gauge group which transforms the gauge field along a gauge orbit, and is by nature non-local. Its 
role in constructing a gauge fixing which is globally consistent and devoid of Gribov ambiguities, 
[27], has been discussed in, for instance, [281 ISHl [301 EI]- That non-locality should play a role 
in aiming to describe infrared gluon dynamics should come as no surprise in that asymptotic 
freedom indicates that only at high energies are quarks and gluons effectively free whilst being 
hard to separate at low energies with lower energy interactions needing to be communicated 
over large distances. For [1| the main initial technical hurdle to be overcome was the fact that 
the non-local operator is an infinite coupling constant series in an arbitrary gauge. Hence to 
do a full perturbative analysis and construct a gauge invariant effective potential was initially 
impossible. However, by taking the point of view that such an effective potential exists then it 
seemed sensible to consider it in one gauge. Specifically, the Landau gauge was chosen whence 
the gauge invariant non-local operator truncates to the one local term [3]. Moreover, 

this operator is renormalizable leading to the successful analysis of the operator's condensation 
in the two loop effective potential, [H [25]. Interestingly, the operator in the Landau gauge 
does not possess an independent renormalization since its anomalous dimension is the sum of 
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the gluon and Faddeev- Popov ghost anomalous dimensions, [32l [33l [M! • A similar structure is 
present in the analogous operator in the Curci-Ferrari gauge, [35l[3B], and the maximal abelian 
gauge, [371 EH]. 

In light of these observations one would still wish to handle the gauge invariant non-local 
operator itself within the context of quantum field theory with the ultimate aim of constructing a 
gauge independent effective potential in order to study the condensation, similar to [4l[25]. As a 

first stage in such an exercise, the main purpose of this article is to renormalize ^uf J d'^x {A'^^)^ 
at one loop in an arbitrary covariant gauge and verify its independence of the gauge parameter 
as well as the equivalence of its anomalous dimension to the Landau gauge expression. From the 
nature of the operator this is not a trivial task as one has to renormalize a non-local operator 
in a gluon 2-point function. Interestingly it transpires that a part of the calculation has already 
been performed but in a different context, \39\ I40|. which we naturally adapt for the present 
work. Before we discuss this it is necessary to indicate the relevant background to the problem 
and indicate the pitfalls of an initial naive approach. First, we set notation for building the 
gauge invariant non-local operator by defining 



O = 



1 



mm 

{u} 



d'x [ai' 



(1) 



Given a gauge field = A'^T"', where T'^ are the colour group generators, it is transported 
along a gauge orbit by the (global) gauge transformation 



UAM^ 



(2) 



where g is the coupling constant. By construction, the field A^ is gauge invariant and hence it 
is trivial to see that the operator O is gauge invariant. However, the explicit form of can be 
determined by setting 

U = e*'^"^" (3) 
where (j)"" can be deduced order by order in perturbation theory. For instance, [301 139] . 
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where are the colour group structure constants, from which it follows that 



O 



1 1,'. 
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+ Oig' 



A? 



(5) 



In QED the operator truncates and represents the square of the transverse component of the 
gauge field. In the non-abelian case the operator involves an infinite number of terms and is 
assumed to converge. 

Naively one would expect to be able to renormalize this version of O order by order in 
perturbation theory since at any order there are only a finite number of terms. It transpires 
that this is not possible without generating several difficult technical problems. First, at one 
loop one would have to extract the divergent part of the topologies represented in figures 1 
and 2 where the encircled cross denotes the location of the operator insertion. Ordinarily to 
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Figure 1: Two leg operator insertion in gluon 2-point function. 



renormalize the graph of figure 1 one inserts the operator at zero momentum. However, as is 
well known, (see, for example, [H]), one can obtain spurious results since the basis of operators 
into which the Feynman integral decomposes is not closed or complete for this momentum 
configuration. Therefore, one has to have non-zero momentum operator insertion. Nullifying 
one of the remaining external momenta, though, results in additional infrared divergences which 
need to be handled. To circumvent these difficulties one approach is to introduce a spurious 
infrared regularizing mass which allows for the nullification of external momenta. Whilst this 
is in principle possible for local operators, for the particular operator we are concerned with it 
has an inherent non-locality which could lead to further difficulties at higher loops. 



Figure 2: Three leg operator insertion in gluon 2-point function. 

Rather than trying to handle these technical issues, it seems more appropriate to regard O 
in a different way and use current results to extract its anomalous dimension. In [30j the gauge 
invariant operator O was rewritten as the sum of gauge invariant operators which can be treated 
individually in the renormalization procedure. Specifically, writing the summation of ([5]) in this 
non-perturbative way, |30j, we have 



where 
and, [30] , 




with the covariant derivative, D^, and field strength, G^j^, given by 

D,Al = d,Al - gr'^A'^Al , G% = d,Al - d.A^ - gr'^A'^A^ . (9) 

Indeed one could regard this as an expansion in operators with n legs where n is determined 
from the lowest number gluon legs in each sub-operator On- 
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The first term ([7]) of ([6]) has already been studied in depth in [39j|40j where it was considered 
as a potential alternative mass operator for the gluon which was gauge invariant though non- 
local. Its anomalous dimension has been computed to two loops in the MS scheme in an arbitrary 
linear covariant gauge and is given by, [391 HO] , 

1 \ abed \ abed „ 

lo,{a) = - -[llCA-^TFNf]a - — + f-^^ f'^^ \-''^- 



3 ' 128Na ' ' 8Na 



+ 



a' + O(a^) (10) 



in terms of a = (^^/(IGvr^) where 

tr(r"r'') = Tpd"'' , T^T" = Cpl , f^'^^f'"^ = CaS^'' (11) 

and Na is the dimension of the adjoint representation. Clearly, the final expression is indepen- 
dent of the usual gauge parameter. To determine this result, the approach was to first localize 
O2 by introducing a set of additional localizing fields and associated (anti-commuting) ghost 
fields in such a way as to produce a renormalizable operator, |39l I40j . From the algebraic renor- 
malizability analysis to ensure a multiplicatively renormalizable localization additional quartic 
interactions between all the localizing fields with couplings, X"-^^'^^ are required. These appear 
for the first time at two loops in the operator anomalous dimension, (llOp . Therefore, in the 
context of determining the anomalous dimension of O the contribution from the first term of dH) 
is known. As we are only concerned with one loop, the piece involving the quartic couplings will 
not become relevant before any two loop renormalization. Hence, to complete the calculation of 
the anomalous dimension of O, 70(0), all that is required is the piece deriving from ([8|) whose 
contributions will be deduced from the graphs of figure 2. These give 703(0) whence 70(0) 
emerges at one loop from the sum of 702(0) ^"^^ 10z[o)- 

At one loop this is actually a simple calculation primarily as a result of the topology of the 
two graphs of figure 2. Since the operator connects to an external leg, the problem of whether 
a zero or non-zero momentum actually flows through the operator insertion does not arise. The 
net flow through the combination of external leg and connecting operator is non-zero which can 
be distributed across both with neither being zero. In other words none of the earlier external 
momenta nullification complications arising in the two leg insertion of figure 1 will arise for the 
graphs of figure 2. Moreover, the non-locality resident in the operator does not lead to any 
additional difficulties and the Feynman rule is simple to derive. This is due to the fact that 
the operator l/D"^ can be replaced by 1/9^ at one loop and this only acts on one field for the 
leading leg term. Therefore, we have computed the two graphs of figure 2 with inserted in 
an arbitrary covariant gauge and found that the sum of the contributions to 70 (a) from both 
graphs is 

703(a) = -^CAa + 0{a^) (12) 

independent of the gauge parameter. For this particular calculation, we used the Mincer 
algorithm for the evaluation of massless 2-point Feynman diagrams, [42^ I43j . written in the 
symbolic manipulation language Form, [31], where the electronic representation of the graphs 
were generated via the Qgraf package, [l5], before being converted to Form input notation. 
Moreover, we use dimensional regularization in d = 4 — 2e dimensions and absorb the divergences 
into the renormalization constants using the MS scheme. The validity of using the Mincer 
algorithm for massless propagators follows because the Feynman integrals are infrared safe and 
no infrared regularizing mass needs to be introduced. The respective numerical contributions to 
(jl2p from the two operators of ([8]) are 33C^/4 and — 15C^/2, which are each independent of 
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the gauge parameter since each operator is itself gauge invariant. Thus adding all contributions, 
we find 

7o(«) = - -^[35CA-l6TFNf]a + ©(a^) (13) 

where the independence of the gauge parameter is a trivial consequence. Moreover, the result is 
in agreement with the anomalous dimension of computed in the Landau gauge, [301 [311 

I32j . Additionally the result appears to be consistent with the expectation that the anomalous 
dimension of a local gauge invariant operator is gauge independent, even though the operator 
itself is non-local. 

We conclude by noting that it would appear that one can renormalize the non-local operator 
O order by order in perturbation theory using Zwanziger's non-perturbative decomposition into 
gauge invariant non-local operators. However, a more stringent check of this would require 
a full two loop calculation which is a non-trivial exercise. Although the contribution to the 
overall anomalous dimension from the localized version of the two leg part of the operator is 
already available, [39\ HOj . much of the groundwork has to be developed for the three and four 
leg operator insertions, such as the Feynman rules. Moreover, the four leg gauge invariant 
operator has yet to be constructed. In addition, in order to apply the MiNCER algorithm it is 
not inconceivable that both O3 and ©4 would need to be localized first by extending the original 
algebraic renormalization analysis of O2, [391 I40j . Nevertheless, such a calculation would be 
interesting since it could validate the observation that the full anomalous dimension of O is in 
fact given by the Landau gauge value which in turn is determined from the sum of the gluon 
and Faddeev-Popov ghost anomalous dimensions in that specific gauge. 
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